function [acceleration_fft] = my_fft(acceleration, filename)

acceleration_fft = fft(acceleration);
fs = 50;% sample frequency = 1 / 0.02 = 50 Hz
N = length(acceleration);% number of sampled points
n = 0 : N-1;
f = n * fs / N;% frequency
plot(f, abs(acceleration_fft));
% saveas(gcf, filename, 'bmp');

end

